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Q ' Lee-Suzuki similarity transformations and Krencigfowa-Kuo folded diagrams 

2^ ■ are two common methods used to derive energy independent model space 
effective interactions for nuclear many-body systems. We demonstrate that 

"^ , these two methods are equivalent to a Renormalization Group (RG) analysis of 
a well-studied problem in quantum mechanics. The effective low- momentum 

CN ' potentials Veff obtained from model space methods are shown to obey the 

r~^ ■ same scaling equation for Vejf that RG arguments predict. This indicates that 

r^ ■ model space methods might be of interest to those studying low-energy nuclear 

^iT . physics using Effective Field Theories (EFT). We find the new result that all 

O ! of the different energy independent model space techniques yield a unique 
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low-momentum 14// when applied to the toy model under consideration. 
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I. INTRODUCTION 

There has been much effort over the past decade to describe low-energy nuclear phenom- 
ena using the techniques of Effective Field Theory (EFT) |[1|-0J^^. The main goal of an EFT 
description of these phenomena is not to out-perform the traditional methods at fitting the 
experimental data. Rather it is to provide a model independent effective low energy theory 
whose form is constrained by the symmetries of QCD, as well as to eliminate the uncon- 
trolled and unjustified approximations that are made in the traditional approaches. One of 
the more appealing features of an EFT treatment is the ability to reliably estimate the errors 
in calculations. This allows one to consistently calculate observables to any chosen level of 
accuracy in principle. The success of various low energy effective theories (Landau Fermi 
Liquid theory, the Fermi theory of Beta decay, the Standard Model, etc) results from the 
fundamental tenet of EFT: low energy observables are insensitive to the details of the high 
energy dynamics. Even if the high energy dynamics are unknown, one can still mimic their 
effects on the low energy physics with local, model-independent operators that are consis- 
tent with the low-energy symmetries. The corresponding coupling constants then implicitly 
contain the effects of the integrated out high energy degrees so that the low-energy physics 
of the full theory is preserved. EFT's are founded upon the ideas first put forward by Wilson 
in his formulation of the "modern" renormalization group in the study of phase transitions. 
The early review article of Wilson and Kogut ^ and the more recent papers of Lepage 
and Polchinski |^,^] illustrate these ideas with pedagological examples and a minimum of 
formalism. 

Shell-model theorists have long employed the notion of effective interactions that are 
defined only within a truncated model space, but that implicitly contain the effects of 
the states that are being thrown away (i.e.- integrated out) so that certain low energy 



observables are preserved from the full-theory p!T|-[T6|]. The qualitative similarities between 
the modern RG techniques and model space methods for deriving effective interactions 
within a truncated model space are obvious. It would be interesting to study the quantitative 
similarities between the two approaches along the lines of the recent work of Haxton and 
Song [^], who seek to employ RG concepts to formulate a "rigorous" shell model free of 
ad-hoc assumptions with the ability to reliably estimate errors. The scope of the current 
paper is a bit more modest and is similar in spirit to the earlier work of Fields et. al. ||17[|, 



where they studied the effective interaction for a simple model problem by writing down a 
scaling equation for the "generalized G- matrix" as the UV regulator is varied and eventually 
sent to infinity. Our approach is more in line with the Wilsonian view of renormalization, 
as we keep the UV regulator finite throughout and instead study how K// scales as we vary 
the boundary of the model space. 

The main purpose of this paper is to study a separable model problem that allows closed 
form solutions of the Lee-Suzuki (LS) and the Krenciglowa-Kuo (KK) iteration schemes. 
We find that K// derived from both schemes obeys the same scaling equation that one 
would obtain from a RG analysis of the model problem |^ . Moreover, we find the new and 
interesting result that all of the various model space techniques corresponding to different 
methods of solving the decoupling equation (LS and KK are two specific schemes) give 
the same low- momentum 14//. In other words, the low- momentum Vejf is unique and 
independent of the particular model space scheme one employs for the toy model we consider. 



This is a suprising result as it is known that when one considers bound state (e.g- shell model) 
problems the different model space techniques often converge to different (i.e.- non unique) 

K//'s inHTel. 

II. MODEL SPACE TECHNIQUES 

We turn now to a quick review of the methods one uses to derive energy independent 
model space effective interactions. Andreozzi has shown that both the Lee-Suzuki (LS) and 
Kuo-Krenciglowa (KK) methods can be recast in the form of a general similarity transfor- 



mation [0 . In this formalism we perform a similarity transformation to obtain an effective 
Hamiltonian that acts only within the model space, but that preserves low-energy physics 
(spectra, scattering amplitudes, etc) from the full problem. If we denote the projection 
operator onto the model space as P and its complement as Q (using the eigenkets of the 
unperturbed Hamiltonian as our basis) we have 

S^^HSS-^ I ^) = ES-^ I ^) (1) 

We can write the above equation in an obvious way by defining H^fj = S~^HS and 
I x) = S'^ I \E'). A convenient form of the similarity transformation is iS = 1 + cj, where 
u is called the wave operator and is defined to satisfy to = QuP. It is easy to see that 
S^^ = 1 —uj and that the P-space projections of the transformed and original states are the 
same, P | x) = -P I ^)- We want to exploit the invariance of eigenvalues under similarity 
transformations to preserve a subset of the exact eigenvalues, but with H^ff acting only 
within the model space. In other words we want i^e// = PHejfP and \ x) = P \'^) so that 
the excluded Q-space completely decouples from the problem. A necessary and sufficient 
condition for this is QH^ffP = 0. Explicitly writing this out gives the so-called decoupling 
equation for uj 

uPHQuj + ujPHP - QHQu - QHP = (2) 

The decoupling equation is a non-linear operator equation and different methods of solu- 
tion can give different answers. Andreozzi has recently shown that two iteration methods of 
solving it are equivalent to the Lee-Suzuki method and the Krenciglowa-Kuo folded diagram 
theory |T5[. The differences between the LS and the KK methods are most apparent when 



one is dealing with bound state problems and the basis states are discrete (e.g.- harmonic 
oscillator orbits). If the model space is d- dimensional, then the LS H^jf reproduces the 
lowest d eigenvalues of the full problem and the corresponding P-space projections of the 
exact eigenstates. Conversely, the KK H^ff reproduces the d eigenvalues of the full problem 
whose eigenstates have the largest P-space component (i.e.- maximum overlap). The great 
benefit of Andreozzi's work is that it contains solutions of the decoupling equation that are 
formally equivalent to the LS and KK schemes, but his methods are non-perturbative in 
nature and hence eliminate the need to calculate irreducible vertex functions (i.e.-Q-boxes) 
and perform G-matrix resummations. Andreozzi's methods can furnish analytic results for 
simple toy models, and yet they are extremely robust and easy to implement numerically 
for realistic bare potentials (Paris, Bonn-A, V-18, etc..) []IU|. In deference to Andreozzi's 



simplification of KK and LS calculations, we hereafter refer to them as the Andreozzi- 
Krenciglowa-Kuo (AKK) and the Andreozzi-Lee-Suzuki (ALS) schemes. We simply quote 
the relevant equations and refer the interested reader to Andreozzi's paper for details. 

THE ALS METHOD 

Writing uj = J2'^=o-^n and an = J2m=o-^m, the ALS equations are 

Xn = —, -X„_ip(o-„_i) n = 1, 2, ... (4) 

q\(yn-\) 

p{an-i) = PHP + PHQan-i (5) 

g(a„_i) = QHQ - an-iPHQ (6) 

Note that p{(Tn-i) and q{(Jn~i) are the P-space and Q-space effective Hamiltonians at each 
step in the iteration, and the iteration converges when o"„ ^ (Jn-i- 

THE AKK METHOD 

Using the same notation as the ALS scheme, the X^ in the AKK scheme are given by the 
following equations, 

XqPHP - QHQXq - QHP = (7) 

XnP{an-l)-QHQXn + (Tn-lPHQXn-l = n=l,2,... (8) 



HI. SCHEMATIC MODEL WITH SEPARABLE POTENTIAL 

Both the ALS and AKK methods are easy to implement numerically for realistic models 
of the bare Vnn such as the Paris and Bonn-A potentials. However, it would be nice to 
construct a toy model for which the ALS and AKK equations yield analytic solutions. For 
this purpose, we consider the following schematic underlying (or full space) Hamiltonian 

H = Ho + g\v){v\ (9) 

where 

{i\ H \j) = eiSij + grii7]j (10) 



We will use the following formula many times in the analysis that follows. 

5ij gr]ir]j 1 



^\H-'\j) = ^- """ (11) 



where 



F = T.- (12) 

I ^i 

Note that we don't have to worry about the above expressions becoming singular (e^ = 0), 
as we are always inverting in Q-space (where the Cj 7^ 0). 



IV. ALS AND AKK ANALYSIS OF THE SCHEMATIC MODEL 

For the remainder of this paper, all state labels refer to relative momenta, and it is 
assumed we are working with a particular partial wave (labels are suppressed). If we assume 
that the underlying potential V{k' , /c) ~ for k,k' > A (either by a form factor or a 
sharp cut-off), then the truncated model space will consist of all states with k < -, where 
s > 1. The excluded Q-space will then consist of all states with momenta lying in the shell 
— < k < A. For a finite shrinking or decimation (e.g. s = 2), we can use the formulas of the 
previous two sections to calculate the Xn of both the ALS and AKK schemes analytically. 
Unfortunately the expressions quickly become unwieldy forn > 1 , and it is not at all obvious 
how one can sum all of the X„ to obtain an analytic expression for the wave operator (and 
hence for K//). The trick that gets us off the hook is to consider infinitesimal decimations 
where - = A — SA. In this way we are able to calculate the wave operator analytically and 
obtain a scaling equation for K;//- Using the definition l^// = -f^e// — -^0, cavalierly ignoring 
factors of 27r, and working in units where e^ = p^, we have 

A 

Veff{k',k) = V{k',k)+ I V{k',q)uj{q,k)q'^dq k,k'<A-6A (13) 

A-5A 

For (5A ~> 0, this becomes a flow equation for Ve/f 

^V^ffik', k) = -A'V^ffik', A)a.(A, k) (14) 

Since we take the limit of 5A — > 0, we can safely ignore contributions to u that are of 0{SA) 
and higher. With this in mind, both the ALS and the AKK equations simplify considerably. 

ALS FOR INFINITESIMAL DECIMATIONS 

A 

Recalling that each intermediate projection operator Q = J q^ \ q){q\ dq ~ 0{5A), 

A-5A 

we can replace all of the p(a"„_i) = PHP + PHQan-i by PHP in the ALS equations. It 
can be shown that at each step in the ALS iteration scheme we can write 



{q I g(a„_i) | q) = €q5q'q + gr]q,r]q (15) 

Where rfqi = r]q/{l + constants) . Consequently, we can use equations 11 and 12 to write 



{q' I r^(^n-i) I g) = ^ - ^ + 0{6A) (16) 

Naively, one might think that the X„ are all infinitesimals of 0{6A) since the ALS equations 
all have an intermediate Q projection operator. This is wrong. One must remember that by 
equation 16 , the operators (QHQ)^^ and g~^(o-„_i) consist of a diagonal part proportional 
to a (5-function plus an off-diagonal part. It is clear that when an intermediate infinitesimal 
projection operator Q acts on the off-diagonal part it results in a term of 0{6A) that can 
be dropped. However, the (5-function portion survives the infinitesimal Q projection and 
gives a finite contribution. Therefore, we can make a second major simplification to the 
ALS equations by replacing [QHQ)~^ and g~^(o"„_i) by {QHqQ)~^. The simplified ALS 
equations can be solved by addition, resulting in the following linear integral equation for u 

" = •^" + 07^'^™^ <"' 

Rearranging terms, explicitly displaying indices, and abbreviating integrals as sums, we 
obtain 

A-5A 

J2 coiA,p')Aip',p) = -eAXo(A,p) (18) 

P' 

Where we have defined the matrix A{p',p) = (cp/ — e\)6prp + grjpirjp. Multiplying from the 
right by A^^ we obtain 



u 



(A,p) = -e^Y.M^^P')^'\p'^P) (19) 



Once again, we make use of equation 11 to obtain an explicit expression for A , viz 



1-1/ / N Vp 9np'% 



A-\P'.P) = -^^ - 7 ^^fP TT— ^ (20) 

Cp - CA (Cp - ca) (ep/ - eAJ 1 + gF 

where we have defined F by the following equation 

A-(5A 2 2 
^ 

Since we eventually take SA — * 0, it is clear that F ^ oo allowing us to ignore the second 
term in the above expression for A~^. Noting that Xo{A,p) = ^^^'^'^p and utilizing the 
simplified expression for A~^ , we obtain our final expressions for the ALS wave operator 
uj{A,p) and the ALS flow equation 

^A,p) = =f^ = 4^£^ (22) 

^ ^ A2 - p2 A2 - p2 V / 



and 



9y f> ^ A^V;//(p^A)V;,,(A,p) 
dX^^ffiP ^P) A^^2 (23) 
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AKK FOR INFINITESIMAL DECIMATIONS 

For infinitesimal model space reductions A ^ A — 6A, we simply state the simplified 
AKK equations as we use the same arguments we used to simplify the ALS equations in the 
previous section. Equation 7 simplifies to 

XoPHP - QHoQXo = QHP (24) 

while equation 8 simplifies to 

XnPHP-QHoQXn = n = l,2,... (25) 

It is easily shown that all of the A„ are zero except for Xq. Hence, the AKK wave operator is 
given by Xq and it is found to be the same as the ALS wave operator. We can generalize 
the above results to cover any method of solving the decoupling equation, by once again 
considering an infinitesimal A — i> A — SA. In this case the non-linear decoupling equation 
simplifies to 

uPHP - QHoQu - QHP = (26) 

It is easily shown that the solution to this equation is the same as the ALS and AKK 
solutions. Hence, when we shrink the boundary of P-space by an infinitesimal amount, all of 
the various methods of solving the decoupling equation give the same equation for d^Veff. 
We shall find that we obtain the same scaling equation if we cut off all loop integrals 
in the Lippman-Schwinger equation and demand the observable scattering amplitude be 
independent of A; this is the RG method described in the next section. We can now invoke 
the semi-group property of the RG and make a much stronger statement about our model 
space techniques. Even for finite A -^ -, the different methods of solving the decoupling 
equation give a unique low-momentum Veff for the toy problem under consideration. To the 
best of our knowledge, this is a new and interesting result. This result brings out yet another 
similarity between model space techniques and the Wilsonian RG. In the application of the 
Wilsonian RG to critical phenomena, there are many different methods of integrating out 
the short distance physics to obtain an effective long-wavelength theory. For example, there 
are infinitely many prescriptions for performing block spin transformations and assigning 
values to the block variables for the Ising model. However, the differences between the 
various coarse graining procedures are eventually washed out and one always converges to 
the same long-wavelength effective theory. 

V. RENORMALIZATION GROUP ANALYSIS 

Since we claim that the model space techniques of nuclear many-body physics give the 
same scaling equation for Ve/f that a renormalization group analysis predicts, it is conve- 
nient to quickly review exactly what we mean by a "renormalization group analysis". In 
relativistic quantum field theory when one calculates physical observables perturbatively as 
a power series in the coupling g, one finds that the coefficients of the higher order powers 
of g are divergent. These infinities are the result of loop integrals that diverge at high 

7 



momentum, or equivalently at short distances. One can isolate these infinities by regulat- 
ing the divergent integrals with a UV cut-off or any other convenient method (dimensional 
regularization, Pauli- Villars, etc.), but then one is faced with the unpleasant result that 
calculated observables depend on the arbitrary regulator mass A that has been introduced 
into the problem, as well as the particular method of regularization. As is well known, the 
way out of this dilemna is to notice that for certain field theories (the so-called "renor- 
malizable" theories) we can add counterterms to the theory that have the same operator 
structure as the terms present in the initial Lagrangian. If we are sufficiently clever, we can 
hide the regulated infinities in the coefficients of the counterterms so that when we calculate 
graphs using the new Feynman rules (that are induced by the counterterms), everything 
is finite. Hence, one can absorb the divergences by redefining the couplings and making 
them A-dependent in such a way that the physical observables are A-independent and fi- 
nite as A — >■ cxo. Thus, a change in the UV cut-off scale is compensated by a simultaneous 
change in the couplings and masses leaving the physics invariant. This is the essence, albeit 
a grotesquely simplified description, of the renormalization group in QFT designed to hide 
the infinities at short distances. But by allowing A ^ oo we are assuming that the dynamics 
encoded in the original Lagrangian holds true at arbitrarily large energies. In the language 
of old fashioned perturbation theory, allowing loop momenta to go up to infinity corresponds 
to summing over intermediate states of arbitrarily large energy and momentum. But this 
presupposes that the intermediate states we are summing are devoid of any " new" physics 
at arbitrarily small distances, which is clearly a bold assumption (supersymmetric partners? 
quark and lepton substructure? strings?) since 1 TeV 7^ 00! About the only thing we can 
say for certain about the high energy intermediate states is that they are highly virtual and 
thus localized in space-time. 

The modern approach to renormalization views A as a physical scale at which unknown 
physics comes into play, and not as an artifact that is to be taken to infinity at the end of the 
calculation. In this way we don't sum over intermediate states we don't fully understand. 
Yet we can model the effects of these excluded states on the low energy physics by recalling 
that their virtual nature means that they are propagated over very small distances and 
times (i.e.- they are localized). We can include their effects on the low energy physics by 
writing down every possible local interaction that is consistent with low energy symmetries 
(model independence!), and we can parameterize our ignorance of the the excluded high 
energy states inside of the tunable couplings that multiply each local operator. Although 
there are an infinite number of local operators consistent with the symmetry requirements, 
they all scale as inverse powers of A allowing one to truncate the effective Lagrangian 
to any chosen level of accuracy. The corresponding couplings can then be obtained from 
a finite number of independent experimental measurements, giving the theory predictive 
power to the chosen level of accuracy. In this way the effective theory is consistent with 
the underlying symmetries, the errors are reliably estimated, different quantities can be 
calculated consistently to the same level of accuracy, and one can infer the presence of " new 
physics" when the effective theory fails. 

We can employ similar ideas in analysing our toy potential model. Here, we can impose 
a UV cut-off A in the loop integrals of the Lippmann-Schwinger equation to reflect our 
ignorance of what is really going on at short distances. We then demand that 1^// depend 
on A in such a way that the low energy observables (i.e.- the scattering amplitudes) are 



independent of A. For our separable toy model we can write the half-onsliell (HOS) T- 
matrix as 



1 1 1 f qW, 



{k' I T{u = k^) I k) gr]k'rjk rjk'r]k ^ k'^ - q^ 



dq (27) 



Since we want the physical quantities to be independent of this UV cut-off A, we allow the 
coupling to "run" with A. By setting d\T~^ = (and noting that all the A-dependence is 
inside g) , we obtain a scaling equation for Vejf 

dX^^ffi^ ' ^> A^^P ^^^^ 

This is same scaling equation for Veff that we derived using model space methods. We also 
note that this is the same equation obtained by Birse et. al. |0, except that their Vg/f 
depends on the energy as well (i.e. Veff{k',k) -^ Veff{k',k; E)). We therefore conclude 
that model space techniques preserve the low energy HOS T-matrix and are equivalent to 
renormalization group methods when applied to the separable model. 



VI. CONCLUSIONS 

We have found that model space methods and renormalization group techniques give the 
same scaling equation for the low- momentum Veff. This was accomplished by considering 
a separable potential model in which the model space techniques yield analytic expressions. 
Moreover, we have found that all methods of solving the decoupling equation give the same 
Veff for the toy model under consideration. To the best of our knowledge, these are new 
and interesting results. The generalization of these results to realistic V/vat's will appear in 
a forthcoming paper. 
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